Quasiparticle relaxation dynamics in superconductors with different gap structures: 
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■ Photoexcited quasiparticle relaxation dynamics are investigated in a YBa2Cu307_i superconduc- 

| tor as a function of doping 8 and temperature T using ultrafast time-resolved optical spectroscopy. 

A model calculation is presented which describes the temperature dependence of the photoinduced 
quasiparticle population n pe , photoinduced transmission AT /T and relaxation time r for three dif- 
ferent superconducting gaps: (i) a temperature-dependent collective gap such that A(T) — > as 
Xj^ • T — > T c , (ii) a temperature-independent gap, which might arise for the case of a superconductor 

with pre-formed pairs and (iii) an anisotropic (e.g. d-wave) gap with nodes. Comparison of the 
theory with data of photoinduced transmission \ AT/T\, reflection [AlZ/TZl and quasiparticle recom- 
bination time r in YBa2Cu307_a over a very wide range of doping (0.1 < 8 < 0.48) is found to 
give good quantitative agreement with a temperature- dependent BCS-like isotropic gap near opti- 
mum doping (8 < 0.1) and a temperature-independent isotropic gap in underdoped YBa2Cu307_£ 
O ' (0.15 < 8 < 0.48). A pure d-wave gap was found to be inconsistent with the data. 
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Oh 1 I. INTRODUCTION 

Ultrafast optical time resolved experiments performed on high-T c superconductors in recent years have shown that 
3 ■ a substantial transient change of the optical transmission or reflection can be induced in these materials by ultrashort 
laser pulse photoexcitation (PE) . What makes these studies particularly interesting is the fact that in optimally doped 
superconductors where T c is a maximum, the amplitude of the observed photoinduced transmission signal, \AT/T\ , 
(or reflection \A1Z/1Z\) appears to increase dramatically below T C EJ~Q. Furthermore, recently it was found that in 
underdoped YBa2Cu307_,5 (6 > 0.15), the increase in amplitude could not be correlated with T c , but rather with the 
' so-called "pseudogap" temperature T* which increases with increasing <5tl These observations together suggest rather 
1-1 \ strongly that the photoinduced effects are related to the opening of a gap (or pseudogap) in the density of states close 

• to the Fermi energy. For optimally doped cuprate superconductors - on which most experiments were done so far - 
there has been substantial agreement in the literature regardhMydie experimental data. However, various groups have 

• proposed very different explanations for the observed effectsQ&QtJand their origin have remained controversial]. 

£f) ' In this paper we present a calculation of the photoinduced optical response in a superconductor for weak excitation 
and compare our results with systematic experimental data on YBa2Cu307_a as a function of S and T . We derive 
ON ■ expressions for the temperature dependence of the photoexcited quasiparticle density n pe and photoinduced trans- 
mission amplitude \ AT/T\ for different possible gaps (or pseudogaps) which might be applicable to high-T c cuprates. 
The simplest case considered is that of a gap, whose temperature dependence is mean-field like and is approximated 
using a BCS function. The second case considered is for a temperature- independent "pseudogap" which is relevant 
for a superconductor where condensation of pre-formed pairs into a superconducting state occurs at T c but pairing 
takes place above this temperature. In this case the gap is not due to a collective effect, but essentially represents 
the pair binding energy Eb and is - to first approximation - temperature independent. Finally we discuss the case 
of an anisotropic gap with nodes, such as one might expect for a d-wave superconductor. In section V we present 
a calculation of the temperature dependence of the quasiparticle relaxation rate r. Comparisons of the theoretical 
predictions with experimental data are made in each case and a detailed discussion is given in section VI. 

> 
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The time-resolved experiments were performed using 100 fs, 800 nm pulses from a Tksapphire laser using single 
high frequency modulation of the pump at 200 kHz. The YBa2Cu307_,5 thin film samples were grown on MgO or 
SrTiC>3 substrates and annealed in oxygen to obtain different O concentrations S. The T c was measured by measuring 
the AC susceptibility in each case. The transition widths, defined as the temperature where \ drops to 90% of full 
diamagnetism, were typically 1-2 K for the S < 0.15 and 4-7 K for S > 0.15. The oxygen concentration was determined 
from T c , which brings in some uncertainty near optimum doping, since the optimum T c for YBCO is 90 K on SrTiC>3 
substrates and 89 K on MgO substrates, which is 2-3 K lower than in single crystals. The typical film thicknesses were 
100-120 nm and all experiments were performed in transmission through the sample. The typical energy of the laser 
pulse is 0.2 nJ, and the laser spot size typically 100/im in diameter. In estimating the carrier density, the absorption 
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length was assumed to be I ~ 0.1 /im. The CW temperature rise of the superconductor film due to laser heating 
was calibrated using a 30/xm YBCO superconducting microbridge from the same batch and on the same substrate as 
in the optical experiments by measuring the change in resistivity with a 4-probe contact method. This temperature 
offset was then taken into account in the plots shown in Figures 1-6. The error in the sample temperature is thus 
reduced to less than ±2 K. 

The photoinduced transmission AT /T through two different films with T c = 90 and 53 K respectively is shown in 
Fig. 1 at different temperatures above and below T c . There are two decay components in each case, a fast component 
with r ~ 0.3 — 2ps, .and a longer-lived component with tl > 10 ns, which appears as a nearly flat background already 
discussed elsewhereo. The latter will not be analysed here, except that it will be taken into account in the fitting 
procedure. The time-evolution of the traces was fitted (shown by the solid lines in Fig.l) using a model with a 
single exponential decay and a Gaussian temporal profile pump pulse, from which the amplitude of the photoinduced 
transmission, \AT/T\ , and relaxation time of the fast component, r, was determined. The temperature dependence 
of \AT/T\ and r derived from the fits are plotted in Figs. 2, 4 and 5 and will be discussed together with the theory 
in the following sections. 

III. INITIAL PHOTOEXCITED CARRIER RELAXATION 

The initial phase of PE carrier relaxation after absorption of a pump laser photon proceeds very rapidly. After 
the laser pulse excites an electron-hole pair, these PE carriers thermalize among themselves with a characteristic 
time of r e _ e ~ 2-nE' 1 ' w -h ere E i s the carrier energy measured from the Fermi energy Ep. In this thermalization 

process, quasiparticle avalanche multiplication due to electron-electron collisions takes place as long cLS Tg g IS less 

than the electron-phonon {e-ph) relaxation time T e — p h- Electron-phonon relaxation becomes important when the 

quasiparticle (QP) energy is reduced to E = \J 2 ^tJ h " — 30 — 50meV (assuming E F ~ 0.1 — 0.2eV). r e _ ph has been 

determined experimentally for the case of a Yl^CuaOy-j] as well as Bi2Sr2CaCu208+:r and Bi2Sr2Ca2pU30io+; C i 
from relaxation time fits in time resolved experiments using intense laser pulses. Using Allen's formulate the e-ph 
relaxation time for YBCO has been found to be T e _ p /, = sjjhm — 100fs for initial carrier temperatures T e — Ej/C e 

in the range 3000 k! and r e _ p ^ ~ 60fs for T e ~ 410 KB. Here Ej is the energy density per unit volume deposited by 
the laser pulse, C e is the electronic specific heat, A is a constant which characterizes the electron-phonon interaction 
and (w 2 ) is a mean-square phonon frequency. A has been determined from these experiments to be in the range 
0.9 < r A-.< 1 both in YBa2Cu307_^ and Bi2Sr2CaCu208+a; superconductors. In experiments which we are considering 
hereoQu, the laser intensities are significantly smaller, and the photoexcited electronic temperature T e is only few K 
in excess of the lattice temperature, i.e. T e ~ T, which gives an e-ph relaxation time r e _ p ^ ~ lOfs. Under these 
near-equilibrium conditions where the carrier temperature is near the lattice temperature, this value of r e ^ p h can 
be compared to the relaxation time determined from infrared reflectivity measurements. For energies of the order of 
30 — 50 meV, h~/r ~ 3000 cm _1 (r ~ 16 fs)liil in good agreement with the estimated T e - p y L ~ lOfs. In the absence of 
a gap in the low-energy density of states, this is the timescale of the electron -phonon thermalization. 

In the presence of the gap the situation is strongly modified and a bottleneck in the relaxation occurs after t > 10 fs. 
As a result QPs accumulate near the gap, forming a non-equilibrium distribution. Each photon thus creates 30 ~ 40 
QPs given by £ , pump /2A,where E pump = 1.5eV is the laser pump, photon energy and 2A is the energy gap. Because 
the final relaxation step across the gap is strongly suppressecOO the QPs together with high frequency phonons 
(with fko > 2A) form a near-steady state distribution. The QP recombination dynamics of this system is governed 
by the emission and reabsorption of high frequency phonons. Phonons with fko < 2 A do not participate in the direct 
relaxation of the QPs, since in this case the QP final states would lie in the gap. (A quantitative justification for the 
assumption that phonons play a dominant role in the QP relaxation will be given in section V.) The basis for these 
pump-probe experiments is that this near-equilibrium QP population can be very effectively probed using excited 
state absorption using a suitable second optical probe pulse, giving direct information on QP dynamics and the nature 
of the gap itself. 

Since the probe laser photon energy E pro b e is typically well above the plasma frequency in high-T c superconductors, 
we make the approximation that the transition probability for the probe light is given by the Fermi golden rule, 
with photoinduced quasiparticles as initial states and with final unoccupied states well above the Fermi energy in a 
band which lies approximately at Eq ~ E pro i, e . The amplitude of the photoinduced absorption (PA) | AA/.A| is then 
proportional to the photoinduced transmission \ AT /T\ (and in the linear approximation also to \A1Z/1Z\ ) which is in 
turn proportional to the number of photoexcited quasiparticles n pe . The probe signal is thus weighted by the O-Cu 
charge-transfer dipole matrix element and the joint DOS, so \AT/T\ cx —n pe pf \Mij\ where n pe is the photoexcited 
carrier density, pf is the final density of unoccupied states and Mij — (p. A) is the dipole matrix element. Although 
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recent experiments on YBa2Cu30y_5 with different E pro b e show the existence of a resonance for E pro i, e ~ 1.5 eVel, we 
assume here that the adiabatic approximation can be applied, so the photoinduced carriers do not cause any change 
in pf or Mij. In this case the effect of the 1.5 eV resonance is to significantly enhance the sensitivity of the probe. 



IV. TEMPERATURE DEPENDENCE OF THE QUASIPARTICLE DENSITY 



A. Isotropic gap. 

In discussing the theoretical explanation for the observed effects, it is important to make the distinction between 
experiments (i) in which the photoexcited charge carrier density is substantially less thaa.-the normal state carrier 
density, n pe <C n c so the laser makes only a weak perturbation on the superconductoiraSEl and (ii) those where 
n pe ~ n c and photoexcitation is sufficiently strong to close the superconducting gaptra. Whereas for the latter 
case (ii), a theoretical description was proposed by MazirJa, so far there has been no theoretical calculation for the 
photoinduced effects in a superconductor for the case of weak photoexcitationp.(^ML The number of photogenerated 
quasiparticles in low-excitation density experiments is n pe < 3 x 10 _3 /unit cellHoEI. On the other hand, the typical 
quasiparticle concentration in a high-T c superconductor is no = 27V(0)A ~ 0.2 — 0.4/unit cell, where N(0) is the 
density of states at Ep. The number of photoexcited QPs is small compared to the normal state density n pe /no < 10~ 2 , 
so the weak photoexcitation approximation is clearly justified and the photoexcited quasiparticles make only a small 
perturbation of the distribution functions. Assuming that the energy gap is more or less isotropic J^o nodes), we can 
approximate nonequilibrium phonon (n u ) and quasiparticle (/ e ) distribution functions as followgli 

l 



hw q < 2A 



huj q > 2A 



(1) 



fe = 7 — i — r—r (2) 

ex P(fe^) + 1 

where T is the lattice temperature and T is the temperature of quasiparticles and high frequency phonons with 

hw q > 2A. 

We can calculate n pe using Eqs.(l), (2) and by considering the conservation of energy. Assuming A is temperature 
independent and large in comparison to T, the conservation of energy has the following form: 

(n T , ~ n T ) A + (4, - n\) ^J^-^, = £l (3) 

here n T ' , nx is the number of thermally excited quasiparticles per unit cell at T and T respectively, E\ is the energy 
density per unit cell deposited by the incident laser pulse, fi c is phonon frequency cut off and v is the effective number 
of phonon modes per unit cell participating in the relaxation. Taking into account that n pe = [n T i — rir) *C jjt 
and assuming that nx — 2Y(0)fcsTexp(— A/fcsT), the number of photogenerated quasiparticles at temperature T is 
given by: 

„ = £ -lI± (4) 

PC l + lvfe ex P(- A /W 

For the case of a temperature-dependent gap A(T), Eq.(3) will be slightly modified: 



^' - nT) + mmor^T) {n -' - 4) ) (A(T) + kBT/2) = £l (5) 

We approximate ~ 27V(0) v /7rA(T)fc B T/2exp(-A(T)/T)lli with ksT <C A, resulting in a slightly modified 
expression for n pe tJ: 

_ _ £ I /(A(T) + k B T/2) 



2 ' / « T exp(-A(T)A B r) 



N(o)hn c V ttA(t) 
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We note that in Eqs. (4) and (6), the explicit form of n pe ( and hence also |AA/.A| or \AT JT\ ) depends only 
on the ratio fceT/A, showing that the intensity of the photoresponse is a universal function of fc^T/A as long as 
the particular functional form of the temperature dependence A(T) is the same. Another important feature of the 
expressions (4) and (6) is that at T = 0, n pe oc 1/A(0), which allows a determination of A(0), as long as the 
experimental parameters required to calculate £i are recorded sufficiently precisely. 

Since the decay of each phonon involves the creation of two quasiparticles, the number of high frequency phonons is 
proportional to the square of the number of quasiparticles n 2 T> at temperature T . (Taking into account conservation 
of energ , at higher pulse energies when n pe becomes large, we expect a crossover of the photoresponse from a linear 
dependence on laser power, n pe oc / to a square root dependence at very high intensities.) 

To enable comparison of the temperature dependences of the photoinduced transmission given by Eqs. (4) and (6) 
with experiments, we have plotted them as a function of temperature in Figures 2b) and c) respectively. The values for 
the constants for YBa2Cu30y_5 are: v = 10 — 20, iV(0) = 2.5 — 5 eV -1 spin -1 cell -1 , f2 c = 0.1 eV. At low temperatures 
n pe is essentially T-independent in both cases, falling off at higher temperatures. However, from the figure it is clear 
that the high-temperature behaviour is quite different in the two cases. In the case of the T-dependent gap (Fig. 2c)), 
the decrease in \AT/T\ is much more pronounced and starts quite close to T c at T/T c ps 0.8 dropping to zero at T c . 
In the temperature-independent gap case (Fig. 2b)), the photoinduced transmission starts to drop at much lower 
temperatures near T/T* w 0.4, dropping exponentially at high temperatures. Also, a notable prediction of Eq. (6) for 
the T-dependent gap is a small peak near T/T c « 0.7, which is not present in the case of a temperature-independent 
gap formula (4). 

In Fig. 2a) we have plotted the data for the photoinduced transmission amplitude for YBa2Cu30y_5 as a function 
of temperature for three different S. Both the magnitude and the temperature dependence are seen to be strongly 
doping-dependent. For comparison with theory, superimposed on the theoretical curves in Figure 2 b) and c) we have 
plotted the normalised photoinduced xefiectivity and transmission data for YBa2Cu307_<5 with different 5. The data 
for a number of underdoped samplesO in the range 0.15 < S < 0.48 have been scaled onto a common temperature 
scale as suggested by the formula (4) and plotted as a function of T/T* in Figure 2b). T* is defined as the point 
where the amplitude of signal \AT/T\ drops to some fixed value (for example to 5% of the maximum value). The 
data can be seen to fit the theoretical curve remarkably well. 

For near-optimally doped YBa2Cu307_a (5 < 0.1) the data are scaled by T c and are plotted as a function of reduced 
temperature T/T c . However, this time the data are superimposed on the prediction for the T-dependent BCS-like 
gap Eq. (6) in Figure 2c). In this case also, the fit is seen to be good-including the small maximum at T/T c w 0.7. 
This maximum is particulary well observed in the data by Han et al.El and by Mihailovic et al.El, although it is not 
evident in the data of Stevens et al.Q. 

The experimental data thus appear to scale onto one of the two theoretical curves, the near-optimally doped 
YBa2Cu307_5 data (5 < 0.1) agreeing very well with Eq. (6) derived for a temperature-dependent gap, while the 
underdoped data (5 > 0.15) fit expression (4) for a T-independent gap. The fact that the underdoped sample data 
scale onto one universal curve in Figure 2b), while the near-optimally doped data scale onto the curve in Fig. 2c), is 
a consequence of the scaling properties of Eqs. (4) and (6) and is a rather remarkable confirmation of the theoretical 
model by the experiment. From the fits of the underdoped sample data, we obtain gap values of 2A — (5 ± 1)/cbT* 
for underdoped samples and 2A-~ (9 ± l)fcsT c for the near-optimally doped samples, the latter in good agreement 
with other optical experimentsElrE3. 

Further confirmation of the model comes from the predicted linear intensity dependence in Eqs (4) and (6). In 
Figure 3 we plot the T-dependence of \AT/T\ for near optimally doped sample (T c = 89K) for three different pump 
intensities. The data are seen to scale linearly with intensity as shown by the insert where £i is changed over one 
order of magnitude. (Unfortunately, measurements over a larger range of intensities were not possible because of 
steady-state laser heating.) 

B. Anisotropic gap with nodes 

If the gap is anisotropic and contains nodes on the Fermi surface - such as for d-wave or strongly anisotropic 
s-wave pairing - there is no clear gap in the spectrum. Nevertheless the quasiparticle DOS still has a strong energy 
dependence and a separation of low and high energy quasiparticles still exists, albeit is less pronounced than in the 
isotropic case. For an anisotropic gap with nodes we can approximate the quasiparticle DOS as a function of energy 
e as: 

N(e) = N(0)(-£-J ( £ «A a ), (7) 
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where A a is a characteristic energy scale separating the low and high-energy QPs, and the exponent 77 depends on 
the topology of the nodes on the Fermi surface. For a two dimensional Fermi surface (with nodes) 77 = 1, while for a 
3-dimensional case 77 = 2. 

Using equation (7) we can estimate the average energy accumulated by quasiparticles after each laser pulse: 



AE = 2FN{0){k B T' (r >+ 2 '> - fc s T(" +2 ))/A2 = £ I 



f 

Jo 



3.77+1 



exp(x) + 1 



dx 



(8) 



(9) 



F = 1.80 for 77 = 1, and F = 5.68 for 77 = 2, and T and T are the QP temperature and lattice temperatures 
respectively. At low temperatures, the phonon contribution to the total energy is exponentially small (Eq.(8)) and 
can be neglected. In this case, the high and low energy phonons are separated by an energy scale of the order of A a . 
(Note that in general A a scales with the gap, but is not equal to the amplitude of the gap.) 
Solving Eq.(8) with respect to the QP temperature T' we obtain: 



T' 



2FN(0) + B 



1/(77+2) 



The number of thermally excited quasiparticles is determined by the equation: 

7i T = 2GN(0)k B T^ +1 ^ /AJJ, 



(10) 



(11) 



where G is given by the Eq.(9), except that 77 + 1 is replaced by 77, so G = 0.82 for r/ = 1 and G = 1.80 for 77 = 2. 
Combining the two equations (10) and (11), we derive the following equation for the number of photogenerated 
quasiparticles as a function of temperature: 



2G N(0)k B T < - r '+ 1 ^ 



1 



2FN(0)k B T(v+2) 



rj + 2 



1 



(12) 



The expression for the photoinduced carrier density Eq.(12) is plotted in Figure 4 with 77 = 2 (solid line) and the 
same values of N(0) and £j as before. Instead of an exponential fall-off with increasing temperature, the photoexcited 
quasiparticle number decreases rapidly well below T c according to a power law. In the 2-dimensional case with 77 = 1, 
(and the same parameters), the temperature dependence is shown by the dashed line in Figure 4 and differs only at 
the lowest temperatures. Using a T-dependent gap (e.g. of BCS form) in Eq. (12), the curve for n pe is virtually 
indistinguishable from the case where AJJ is temperature- independent. The reason for this is that the effect of the 
T-dependent gap is only important as T — ► T c , but there 77, pe is already small. 



At low temperatures we can neglect the 1 in the round brackets and obtain: 



2GN(0) 



2FN(0) ) 



(13) 



A feature of the d-wave model - which is particularly important for comparison with experiment - is the peculiar 

■7 + 1 

sub-linear behaviour of n pe with laser pump intensity, n pe <~ £j +2 at low temperatures. In the high temperature limit 
on the other hand we can expand the term in brackets to obtain a crossover to linear behaviour, with a temperature- 
dependent slope: 



n pe = 2G^-^- 
k B T 77 



(14) 



It should be pointed out that in the high temperature limit, the phonon contribution should become increasingly 
important and the decrease of the number of photogenerated QP will eventually become exponential, just as in the 
case of an isotropic gap. 

The data for the temperature dependence of the induced transmission for underdoped and optimally doped 
YBa2Cu 3 07_,5 are also plotted on Figure 4, superimposed with the predicted response for the d-wave case. Clearly 
the T-dependence data cannot be described using a pure d-w&ve gap irrespective of dimensionality or the form of 
temperature dependence of A(T) (Eq.(12)). The linear intensity dependence shown in Fig. 3 is consistent with this 
observation, and is also not consistent with a d-w&ve gap, which predicts a sub-linear response. We note that although 
the data apparently exclude the possibility of a pure d-w&ve gap in underdoped and near-optimally doped YBCO, it 
does not entirely rule out a mixed symmetry gap with some <i-wave admixture. 
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V. THE QUASIPARTICLE RELAXATION RATES. 



As discussed in the section III, the relaxation rate of the photoinduced QPs is dominated by the energy transfer 
from high-frequency phonons with Huj > 2 A to phonons with fku < 2 A. To describe the relaxation of nonequilib- 
rium quasiearticles in this case we consider the kinetic equation for phonons, taking into account phonon-phonon 
scattering?!^ 

= Iph-ph{riu} (15) 
where phonon-phonon scattering integral has the form: 

I P h- P h{n w } = 2ir ^ \w quq2 \ 2 {- [(n w + 1) n Ul n W2 - n u (n Ul + 1) (n W2 + 1)] 6(u — ux- uj 2 ) 

+ [(n u + 1) (n Ul + 1) - n u n Ul (n«, 2 + 1)] 5(u) 2 -uj- ui)}, (16) 

w 9i,92 i s the anharmonic coupling constant and is the phonon distribution function. (Conservation of the phonon 
momentum is not relevant because of umklapp scattering.) We neglect the electron- phonon scattering integral, because 
at normal temperatures the relative contribution of quasiparticles to the relaxation of phonons is small due to the 
relatively small number of quasiparticles compared to the phonons N(0)£l c <C v. Electron-phonon collisions can play 
an important role only at temperatures T <^T C but since - as shown in the previous section - n pe is nearly constant at 
low temperature, T cannot be small. The distribution function Eq.(l) reduces to the parts of the collision integral 
Eq.(16) which describe the scattering of high and low frequency phonons separately. Only the part of the collision 
integral that does not conserve the number of phonons with uj < 2A differs from zero. 

To estimate the relaxation time, we multiply Eq. (15) by uj q and sum over q, satisfying the condition < ui q < 
2A(T): 

—gf' = ^2 hw q I ph - ph {n u }. (17) 
If we suppose that the coupling constant is momentum independent w qi _ q2 = w, then: 

r2A 



f f 1 

^ Uqlph-phiriu,} = 27rw 2 J uip (uj) duj J {- 



[(n u + 1) njn^_j - (n u > + 1) + 1)] p \J ) p (uj 

[{n u + 1) [nj + 1) n u+u > - n u n u > (n u+u > + l)] p(w')p(u/ + w)}dJ (18) 

here p (uj) = "ivuj 2 /Vt^ is the phonon density of states in the Debye approximation. We restrict ourselves to a discussion 
of the relatively high temperature limit where the phonon distribution function can be replaced by hsT/Hai. The 
energy accumulated by the phonons with hw < 2A is: 

E< = J n u p (uj) hwduj = y (19) 

In the integral Eq.(18) the first term is equal to 0, because it has only phonon distribution functions with fko < 2 A. 
The first non-zero contribution describes the decay of the phonon with Huj > 2A to two phonons with fku < 2A: 

^§™ 2 ((H§7) 3 ^'-T)(2A)' (20) 

The second contribution describes the inelastic scattering of high frequency phonon with creation of one low fre- 
quency phonon with w < 2A: 

^Mw)^'-^ 1 (21) 
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Eq.(21) shows that the inelastic scattering of high frequency phonons is dominant in the energy relaxation of the 
energy if A <C Hfl c and we can neglect I\ . 

Taking into account Eqs.(17), (19-21) we obtain expressions describing the relaxation of the equilibrated QP-phonon 
temperature: 

dT 1 

m=v h {T - T) (22) 



J_ = 97rv 2 w 2 k B T'A(T) 
Tph (htt c ) 3 

Note that Eq.(22) is rather similar to that derived by AllenE§ for the temperature relaxation in the electron- phonon 
system in the normal metals, except that Eq.(22) describes the phonon energy relaxation modified by the gap in the 
quasiparticle spectrum. 

The phonon-phonon relaxation time (Eq.23) can be expressed in terms of an experimental parameter, namely the 
Raman phonon linewidth r w . Using the Fermi golden rule we can calculate for k B T <C hu>: 



The phonon relaxation rate can thus be expressed in terms of I\j in the following form: 

1 12r w fc s T'A(T) 



Tph 



(25) 



For a temperature dependent gap A, therfiekwation time is expected to show a divergence t p h oc 1/A(T) as is 
indeed observed in optimally doped cupratesEna'tllia. A similar divergence of the relaxation time has been calculated 
previously by Schmidt and SchonEj and TinkhamEJ, albeit for somewhat differently created nonequilibrium situations. 

The typical relaxation timescale T p h given by the formula Eq.(25) is very close to the experimentally observed 
values. From the data on the Raman linewidth of the Ai s -symmetry apical 0(4)_phonon mode in YBCO - which 
has been shown to be particularly anharmonic - T w sw 13 cm -1 and u ~ 400 cm _1 Ell. At T ~ T — T c /2, and using 
A Tc / 2 ~ 200 cm -1 we obtain r ph — 0.8 ps. 

At low temperatures the quasiparticle temperature T is much higher than the lattice temperature T and the 
formula is expected to fail. However, from the experiments (Figure 2) we see from the temperature dependence of 
\AT/T\ that the number of photogenerated quasiparticles is nearly constant at low T. T can therefore be estimated 
from the equation: 

n T / -n T = £j/A(0) (26) 
Taking into account that n T ~ 2N(0)A(0) exp(-A(T)/fc B T) we obtain: 

k B T ~ A(T)/ In |l/(£ / /27V(0)A(0) 2 + exp(- A(T)/k B T))} (27) 

At low T, the nonequilibrium temperature T ~ A(T)/ In {(2A r (0)A(0) 2 /^/} - T c /2, while in the high temperature 
limit, the exponent in the logarithm becomes large and T ~ T. Combining Eq.(25) and Eq.(27) we obtain an 
expression for the QP relaxation time as a function of lattice temperature T and photoexcitation energy Ei which is 
valid for all temperatures < T < T c : 

J_ = 12r^A(T) 2 

T ph ^2 ln{1/(£ - //2A , (0)A(0) 2 +exp( _ A(r)/fcBr))} - I > 

To enable comparison of this formula with experiments, we again use Raman data on high-frequency optical phonon 
linewidths, a value of gap frequency of 200 cm -1 (2A ~ 9k B T c ) from the fits in Figure 2c), the same values of 
jV(0) = 2.2 — 5 eV -1 cell -1 spin -1 as in section III, and a BCS function for A(T). £/ was again calculated using an 
incident laser energy of 0.2 nJ incident on a 100 nm thick film with a spot size of ~ 100/im diameter. We obtain a 
temperature dependence of the relaxation ti»e r p h as shown in Figure 5. A comparison of the theory with the data 
for optimally doped YBa2Cu30y_5 samplcaM in Figure 5 shows good agreement with the theoryc 2 ] with no adjustable 
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parameters, especially near T c . Some experimental data in the literatureH'B^I show an upturn in relaxation time at 
low temperatures, which is also weakly present in the data in Fig. 5. Although Eq.(28) is derived for temperatures 
close to T c , an upturn at low temperature can still be reproduced by this formula if we reduce the energy per pulse 
Si. (In fact in the limit of Ei -> 0, Eq. (28) gives r ph oc 1/T at low T.) Indeed by setting £i/N(0) in Eq. (28) as 
an adjustable parameter, rather than inserting fixed values as we have done, the low-temperature upturn in the data 
can be reproduced more accurately. 

To end this section we briefly make a quantitative comparison of the relaxation rates for a gapped and gapless 
material. Comparing the relaxation time for the case of a gapped material with relaxation time in a normal metal 
described by Allen's formulate we obtain: 

- XH3uji (29) 



-ph 



4ir(k B T) 2 AT u 



Inserting experimental values for uj and r w we find the anharmonic phonon decay time in a superconductor to be 
approximately two orders of magnitude slower than electron-phonon relaxation in a gapless material as T p h/T e — P h ~ 
60 - 200 for T ~ T c as expected. 

The calculation was performed assuming that the rate limiting step for QP relaxation is anharmonic decay of high- 
energy phonons, while direct electron phonon relaxation was assumed to be small. This is now justified quantitatively 
by a calculation of the direct electron-phonon relaxation rate. The contribution of phonon-electron scattering to 
the relaxation of phonon system can be estimated by deriving the corresponding relaxation rate r^ 1 using the usual 
electron-phonon collision integral (see Appendix): 

i_ = ^v(o)AA(r)r> c 

Tl 2v 1 ' 

We see that the electron-phonon relaxation rate, when compared to the phonon relaxation rate is reduced by the 
factor N{Q)Ml c /v < 1: 



T ph _ Trh 2 N(0)Xuj 2 n c 



n 2AvT u] k B T' 



< 1 (31) 



and can therefore be neglected. However if the number of phonon modes involved in the relaxation is significantly 
smaller, the ratio T p h/ri may approach unity. In this case the electron-phonon collisions can also contribute to the 
phonon relaxation. However the temperature dependence of the relaxation rate 1/ri is very similar to the anharmonic 
relaxation r p h, both being proportional to A(T) and so the two contributions to the relaxation may be difficult to 
separate experimentally. The data shown in Figure 5 are fit to r p h without any fitting parameters, so it appears to 
be reasonable to assume that anharmonic relaxation described by Eqs. (25) and (28) is dominant in YBa2Cu307_^ 
near optimum doping. 



VI. DISCUSSION 



Near optimum doping in YBa2Cu307_«5 (8 < 0.1) we find that the theoretical model gives a good quantitative fit to 
the temperature dependence of the photoinduced transmission (Figure 2c)), provided we use a temperature-dependent 
gap (Eq. (6)). The relatively square shape of the temperature dependence of \AT/T\ , the weak maximum near 
T/T c — 0.7 and the rapid drop in the signal amplitude above this temperature is a consequence of the particular 
temperature dependence of the gap with the property that A(T) — > as T — > T c . In underdoped YBCO on the other 
hand (Fig. 2b)), the temperature dependence is much weaker and can only be reproduced by the calculation using a 
temperature-independent gap (Eq. (4)). It is not possible to reproduce the temperature dependence of the optimally 
doped data using a temperature-independent gap and vice versa, it is not possible to describe the underdoped data 
with a T-dependent BCS-like gapEi 

A similar conclusion regarding the evolution of the gap with doping is reached on the basis of the comparison of 
the calculated temperature dependence of the relaxation time r with the data. The divergence at T c shown in Figure 
5 for samples near optimum doping together with the temperature dependence of \AT/T\ lead to the inescapable 
conclusion that there is a gap in YBa2Cu307_,5 with (S ~ 0.1) which closes rapidly at T c . On the other hand, the 
notable absence of this divergence of r for underdoped YBa2Cu3 07_a (5 > 0.15)u is consistent with the existence of 
a T- independent gap deduced from the temperature dependence of \AT /T\. 

Plotting the temperature T* and T c as a function of 8 in Figure 6, we find that T* drops with increasing dopiri|g-|iii-a 
characteristic fashion, following the so-called "pseudogap" behaviour reported by other experimental techniquesEa~E3. 
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Near optimum doping, the "pseudogap" becomes observable close to T c . Such a situation occurs in BCS supercon- 
ductors where gap formation is a collective effect and Cooper pairs form simultaneously with the formation of a 
phase coherent collective state at T c . However, from Figure 2c) and Figure 3 there is some evidence of remaining 
QP excitations above T c even near S ~ 0.1. The T-independent "pseudogap" thus appears to be still present near 
optimum doping, and the BCS case is not quite realized. 

The underdoped state is distinctly different, with a T-independent normal state pseudogap appearing at T* ^> T c . 
The changes in the photoinduced transmission \AT/T\ at T* in this case are simply a result of an increase population 
of pairs with decreasing temperature. Since the observed changes in \AT/T\ are associated only with pairing and not 
phase coherence (which occurs at T c ), we_see no change in \AT/T\ at T c in agreement with predictions of pre-formed 
pairing models of the underdoped stateEj'Ej. The pseudogap in this case signifies the pair binding energy, and its 
value is determined by the microscopic pairing mechanism. It is to first approximation temperature independent and 
decreases with an inverse law with increasing carper concentration, which is consistent with increased semi-classical 
(Debye-Hiickel) screening between charge carriersEa. 

Since the temperature dependence of \AT/T\ is actually exponential at high temperatures (Eq.4), it is more 
appropriate to discuss the doping dependence of the magnitude of the energy gap A rather than of T* . This avoids 
the somewhat arbitrary criterion for determining T* of an asymptotic function. Such a plot of A versus doping S 
is shown in Figure 7. As predicted on the basis of screening arguments above, A appears to follow an inverse law 
A oc l/i where x = 0.6 — S is proportional to the carrier concentration. 

The doping dependence of A can be independently confirmed from the same data set by plotting the low- 
temperature value of \AT/T\ as a function of 6. As can be seen from the formulae (4,6), the value of \AT/T\ 
at T = is proportional to £//A, from which the zero-temperature gap A^=o can be determined. A scaled plot of 
1/ \AT/T\ t^o as a function of doping 5 in Figure 7 independently confirms the inverse relation A oc ^JlQ,.6 — S). We 
note that the observation of an inverse law is in agreement with recent experiments on I^-aSr^CuC^Ei^Eil suggesting 
that it may be a universal feature of the cuprates. 

A somewhat surprising feature of the data is that they are described so well by an isotropic gap over the whole 
range of doping. A d-wave gap is apparently not consistent with either the temperature dependence of \AT/T\ or the 
linear intensity dependence of \AT/T\ . It is also not consistent with the observation of a single exponential decay of 
AT /T with time, since in the rf-wave case we would expect a distribution of rs, and certainly no single divergence 
of t at T c should be observed in the d-wave case. A possible explanation for the absence of a d-wave signature is 
that in YBCX) the gap in the bulk of the sample is more or less isotropic in contrast to the surface, where it may be 
more d-likeEa. Another possibility for the apparent discrepancy between the present optical data and, for example, 
Raman experimentsuJ are the different weights that the matrix elements have in the probing process. In the present 
experiments the dipole transition matrix elements average over the whole Fermi surface (FS) weighted by areas of 
large electron momentum. Anisotropy of the FS is not emphasized in this case. This is different than in Raman 
measurements, where the electronic scattering intensity is zero for an isotropic parabolic band, and Raman intensity 
comes only from the anisotropic regions of the FS. 

We conclude by noting that the comparison of the theoretically predicted temperature dependence of the photoin- 
duced absorption signal amplitude and relaxation time r with the experimental results on underdoped and optimally 
doped YBa2Cu307_a gives good quantitative agreement. The evolution of the gap from a temperature-independent 
A to a temperature-dependent collective gap A(T) is a particularly striking feature of the data as is the apparent 
inverse dependence of the gap magnitude with doping. Although the model calculation presented in this paper was 
performed with the superconducting cuprates in mind, it is quite general and can, for example, also be applied to a 
case of photoexcitation experiments on CDW systems or narrow gap semiconductors. It could also be extended for 
other non-equilibrium situations as might occur in non-equilibrium superconducting high-T c devices for example. 
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VIII. APPENDIX 

To estimate the electron-phonon relaxation rate we take into account the electron-phonon collision integra 







Iph-eM = -4vr^ \M kk ,\ 2 [f k (l - f k ,)n q - f k ,(l - f k ) (n q + 1)] 6(e k - e k> + w q ) 

k 

M kk ' is electron-phonon coupling constant and n u , fk are phonon and electron distribution functions. 

Note that we neglect coherence factors in the electron-phonon collision integral. This is a reasonable approximation 
for the BCS case if the characteristic energy of quasiparticles is large in comparison with the gap T > A. For 
underdoped case also, the normal state gap is temperature independent and hepce coherence factors are equal to 
It is easy to show that the collision integrals satisfy the conservation energy lawEJ: 

k q 



To calculate electron-phonon collision integral we define "electron-phonon spectral function" a F : 
a 2 F(e,e',n) = ^ £ |* 5(u> q - Sl)5(e k - e)6(e k , - e) 

k,k' 



Following AllenEj, we suppose that a 2 F(e,e , fi) ~ a 2 F(ep, ep , SI) = a 2 F(fl). For the sake of simplicity we make 
le Debye approximation for the electron-phonon spectral function a 2 F(Q) — 

<>2A poo 

Uqle-phinu,} = -27riV(0) / n a 2 F(n)dn {(/( e )_/( e + n))nn- 

,K<2A) J ° 7A 

/(e + n)(l-/(e))}de (Al) 
The integral Eq.(Al) can be easily evaluated: 



4 



h 2 n 2 

Formula (A2) leads to the Eq.(30) for electron-phonon relaxation rate. 



z s ,i«^(r-r, (A 2 , 
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Figure 1. The normalized photoinduced transmission AT /T as a function of time delay for two different samples 
with a) T c — 53K and b) T c — 90K, each at three different temperatures. The time-evolution of the traces was fitted 
(solid lines) using a model with a single exponential decay and a Gaussian temporal profile pump pulse. 

Figure 2. a) The photoinduced transmission amplitude \AT/T\ in YBa2Cu30y_5 as a function of temperature for 
three different 8 : for 8 — 0.48 (open squares), 8 = 0.18 (solid triangles) and 8 ~ 0.1 (open circles), b) The solid 
line is a plot of the photoinduced transmission amplitude as a function of temperature for a temperature-independent 
gap using Expression (4) as a function of T/T* . The data points are for underdoped samples with 8 = 0.18 with 
T c = 77 K (solid triangles), 8 = 0.32 with T c = 62 K (open triangles), 8 = 0.44 with T c = 53 K (solid circles) and 
6 = 0.48 with T c = 48 K (open squares) respectively, c) The solid line is the calculated photoinduced transmission 
amplitude from Eq. (6) for a temperature-dependent gap as a function of normalized temperature T/T c and data for 
the near-optimally doped samples 6 ~ 0.1 (T c = 90K) from Mihailovic et alH-Zopen circles), Stevens et alia (open 
squares) and Han et alB (solid triangles). Note: the original datapoints of Refill were shifted by 7K in temperature. 
The solid lines in a) are fits with Eq. (4) for 6 = 0.48 and 8 = 0.18, and Eq. (6) for 5 ~ 0.1. 

Figure 3. The T-dependence of \AT/T\ for near optimally doped sample (T c = 89K) for three different laser 
intensities. The data has been normalized with respect to the energy density per unit volume deposited by the laser 
pulse, £j. The insert shows low temperature \ AT/T\ vs. £j, where \ AT/T\ scales linearly with £j and £q = 8 x 10~ 5 
eV cell _1 spin _1 . 

Figure 4. A plot of the temperature dependence of the photoinduced transmission amplitude ( i.e. PE carrier 
density) for a pure <i-wave gap (Eq. (12)). The solid line is for a 3D case, while the dashed line is for a 2D case using 
the same parameters. (The ordinate scale is the same for both curves.) The data (symbols) exhibit qualitatively 
different behaviour than predicted by the ci-wave model, irrespective of doping 8. The data were scaled in the same 
way as in Fig. 2. The open squares represent data for 8 > 0.15 while the solid squares are for 8 < 0.1. 

Figure 5. The relaxation time r as a function of temperature for a temperature-jdependent gap (Eq. (28)). The 
data are for the near-optimally doped samples 8 ~ Orl (T c = 90K) from Han et al.El (full circles) and Mihailovic et 
al.El (open circles). Note: the data points of Han et altU were shifted by 7K in temperature. 

Figure 6. T c (open squares) and T* (solid squares) as a function of i5. T* is defined as the point where the 
amplitude of signal \AT/T\ drops to 5% of maximum of the T-dependent signal (the T-independent contribution, 
which is assumed to be due to Allen relaxation is subtracted). 

Figure 7. The dependence of the energy gap A on doping 8. The open squares correspond to values of A determined 
from fits to the data using Eqs.(4) and (6), while the solid circles are obtained by plotting the scaled magnitude of the 
inverse of the photoinduced transmission at low temperature, 1/ \AT/T\ T ^ , In the latter case, the data were taken 
in a single experimental run, with carefully controlled laser operating conditions to ensure that A can be compared 
quantitatively. The dashed line is a guide to the eye showing Aal/i bahaviour. 
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